NEW CONSTRUCTION OF ALGEBRO-GEOMETRIC SOLUTIONS TO 
THE CAMASSA-HOLM EQUATION AND THEIR NUMERICAL 

EVALUATION 
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Abstract. An independent derivation of solutions to the Camassa-Holm equation in 
terms of multi-dimensional theta functions is presented using an approach based on Fay's 
identities. Reality and smoothness conditions are studied for these solutions from the 
point of view of the topology of the underlying real hyperelliptic surface. The solu- 
tions are studied numerically for concrete examples, also in the limit where the surface 
degenerates to the Riemann sphere, and where solitons and cuspons appear. 



1. Introduction 
The Camassa-Hohn (CH) equation 

(1.1) Ut + SuMa: = Uxxt + 2UxUxx + UUxxx - ^kUx, 



was first found by Fokas and Fuchssteiner |19| with the method of recursion operators and 
shown to be a bi-hamiltonian equation with an infinite number of conserved functionals. 
Camassa and Holm ^12j sliowed that it appeared as a model for unidirectional propagation 
of waves in shallow water, u{x, t) representing the height of the free surface about a flat 
bottom, k being a constant related to the critical shallow water speed. In this context, 
only real-valued solutions are physically meaningful. 

To be able to formulate a Cauchy problem for the CH equation which implies the 



solution of (1.1) for a given function u{x,0), it is convenient to write (1.1) in the non- 
local evolutionary form, 

(1.2) Ut = D~^{-3uux + 2uxUxx + UUxxx - 2kux), 

where the operator D is given by D = 1 — dxx, and where its inverse is defined by giving 
certain boundary conditions. Since the inverse of the operator D stands for an integral 
over the Green's function for given boundary conditions, and thus an integral from some 
base point xq to x, it is not a local operation in x. This non-locality has mathemati- 
cally interesting consequences: the CH equation has traveling wave solutions of the form 
u{x,t) = c exp{ — |x — vt\} {v being the speed, c = const) called peakons that have a 
discontinuous first derivative at the wave peak. Camassa and Holm [12] described the dy- 
namics of the peakons in terms of a finite-dimensional completely integrable Hamiltonian 
system, namely, each peakon solution is associated with a mechanical system of moving 
particles. The class of mechanical systems of this type was further extended by Calogero 
and Frangoise in P1[TU]. Multi-peakon solutions were studied using different approaches 
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in a series of papers [U O [U [TT]. Periodic solutions of the shallow water equation were 
discussed in [29j . 

A further consequence of this non-locality is that solutions of the CH equation in terms 
of multi-dimensional theta functions do not depend explicitly on the physical coordinates. 
Such solutions were first given by Alber and Fedorov in |2j by solving a generalized Jacobi 
inversion problem. In contrast to the well known cases of Korteweg-de Vries (KdV), non- 
linear Schrodinger and sine-Gordon equations, see for instance [7] and references therein, 
complex solutions of the CH equation are not meromorphic functions of {x, t) but have 
several branches. This is due to the presence of an implicit function y{x, t) of the variables 
X and t in the argument of the theta function appearing in the solutions. A monodromy 
effect is thus present in the profile of real-valued solutions such as cusps and peakons. This 
means that even bounded solutions to the CH equation can have discontinuous or infinite 
derivatives in contrast to KdV solutions. Algebro-geometric solutions of the Camassa-Holm 
equation and their properties are studied in [H [21 |3l |22l |23l . 

Our goal in this paper is to give an independent derivation of such solutions based 
on identities between multi-dimensional theta functions, which naturally arise from Fay's 
identity |18) . This identity states that, for any points a,b,c,d on a compact Riemann 
surface of genus g > 0, and for any z G C^, there exist scalars 71,72,73, depending on the 
points o, b, c, d, such that 

(1.3) 71 e(z + o e(z + f^') + 72 e(z + f^^) g(z + x") = 73 e(z) e(z + + 



where G is the multi-dimensional theta function (2.2); here and below we use the notation 
for the Abel map i2.5« between a and b. While the authors in |Tl|2l|3] used generalized 



theta functions and generalized Jacobians (going back to investigations of Clebsch and 



Gordan |T5]), we derive the solutions from the identity (1.3). This fits into the program 
formulated by Mumford |30] that all algebro-geometric solutions to integrable equations 
should be obtained from Fay's identity and suitable degenerations thereof. Historically 
this approach was only able to reproduce solutions already obtained via so-called Baker- 
Akhiezer functions, generalizations of the exponential function to Riemann surfaces. The 
first example of new solutions found via the Fay identity was by one of the authors |27j for 
the multi-component nonlinear Schrodinger equations, see also [28j. Both methods have 
specific advantages: for the Baker- Akhiezer approach, solutions to the associated linear 
system the integrability condition of which is the studied equation have to be constructed 
on a Riemann surface for a given singularity structure. For the Mumford approach the 
non-trivial task is the finding of a suitable degeneration of the Fay identity for the studied 
equation. Once this is done the identification of certain constants in the solutions as 
well as the study of reality ans smoothness conditions is then in general more straight 
forward than in the Baker- Akhiezer approach. We provide here the first example for an 



integrable equation with nonlocal terms in the evolutionary form (1.2) as explained there. 
The spectral data for the theta-functional solutions to CH consist of a hyperelliptic curve 
of the form //^ = rTjfa^(-^ ~ ^j) with three marked points: two of them are interchanged 
under the involution a{X, fi) = (A, —fJ.), and the third is a ramification point (Ajo, 0). 

Our construction of real valued solutions is based on the description of the real and imag- 
inary part of the Jacobian associated to a real hyperelliptic curve (i.e., the branch points 
Xj are real or pairwise conjugate non-real). In this way, one gets purely transcendental 
conditions on the parameters (i.e., without reference to a divisor defined by the solution of 
a Jacobi inversion problem), such that the solutions are real- valued and smooth. It turns 
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out that continuous real valued solutions are either smooth or have an infinite number of 
cusp- type singularities. 

Concrete examples for the resulting solutions are studied numerically by using the code 
for real hyperelliptic surfaces |20|I21|. This code uses so-called spectral methods to compute 
periods on the surfaces. It allows also to study numerically almost degenerate surfaces 
where the branch points collapse pairwise. In this limit, the theta functions break down to 
elementary functions, and the solutions describe solitons or cusps. It is noteworthy that 
the theta-functional solutions thus contain as limiting cases all known solutions to the CH 
equation. The quality of the numerics is ensured by testing the identities between theta 
functions which are used to construct the CH solutions in this paper. In addition, the 
solutions are computed on a grid and are numerically differentiated. These independent 
tests ensure that the shown solutions are correct to much better than plotting accuracy. 

The paper is organized as follows: in section 2 we summarize important facts on Riemann 
surfaces, especially Fay's identities for theta functions and results on real surfaces. In 
section 3 we use Fay's identities to rederive theta-functional solutions to the CH equation, 
and give reality and smoothness conditions. In section 4 we study numerically concrete 
examples, also in almost degenerate situations. We add some concluding remarks in section 
5. 



2. Theta functions and real Riemann surfaces 

In this section we recall basic facts on Riemann surfaces, in particular real surfaces and 
multi-dimensional theta functions defined on them. 

2.1. Theta functions. Let TZg be a compact Riemann surface of genus g > 0. Denote 
by {A, B) := {Ai, . . . , Ag, Bi, . . . , Bg) a canonical homology basis, and by (wi, . . . , ujg) the 
basis of holomorphic differentials normalized via 

(2.1) / ujj = 2m5kj, k,j = l,...,g. 

The matrix IB = ^J^^lvj^ of ;B-periods of the normalized holomorphic differentials lvj, 
j = 1, ... ,5, is symmetric and has a negative definite real part. The theta function with 
(half integer) characteristics 5 = [61,62] is defined by 

(2.2) @m[6]{z)= exp{^(]B(m + (5i),m + 5i) + (m + 5i,z + 2i7r,52)}, 

for any z G C^; here 61,62 G {O, are the vectors of the characteristics 6; (., .) denotes 
the scalar product (u, v) = ^^UiVi for any u, v G C^. The theta function 0[5](z) is even 
if the characteristics 6 is even, i.e., 4 {61, 62) is even, and odd if the characteristics 6 is odd 
i.e., ^{61,62) is odd. An even characteristics is called non-singular if 0[(5](O) 7^ 0, and an 
odd characteristics is called non-singular if the gradient V0[5](O) is non-zero. The theta 
function with characteristics is related to the theta function with zero characteristics (the 
Riemann theta function denoted by 0) as follows 

(2.3) 0[(^](z) = e(z + 2i7r(52 +B(5i) exp { ^(IB^i, (5i) + (z + 2i7r(52, } . 

Denote by A the lattice A = {2i7rN + BM, N, M G Z^} generated by the A and B- 
periods of the normalized holomorphic differentials uoj, j = 1,. . . ,g. The complex torus 
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J := J{TZg) = C^/A is called the Jacobian of the Riemann surface TZg. The theta function 
(2.2) has the following quasi-periodicity property with respect to the lattice A: 

e[5](z + 2i7rN + ]BM) 



(2.4) = G[<5](z)exp{-i(BM,M) - (z,M) + 2i7r(((5i,N) - {52,M))} . 
Denote by 11 the Abel map 11 : TZg i — > J defined by 

(2.5) U{p)= r uj, 

for any p G IZg, where pQ G IZg is the base point of the application, and where uj = 
(cji, . . . , ujgY is the vector of the normalized holomorphic differentials. In the whole paper 
we use the notation = 11(6) — Ii{a). 

Now let ka denote a local parameter near a € TZg and consider the following expansion 
of the normalized holomorphic differentials Wj, j = 1, . . . , g, 

(2.6) UJj{p) = {Va,j + Wa,-j ka{p) + . . .) dka{p), 



for any point p G TZg in a neighborhood of a, where Vaj-, Wa,j G C. Let us denote by Da 
the operator of directional derivative along the vector = (Va^i, ■ ■ ■ , V^,g)*: 

9 

(2.7) DaF{z) = Y,dz^F{z)Va,j, 

i=i 

where F : — > C is an arbitrary function. According to [30], the theta function satisfies 
the following identities derived from Fay's identity (1.3): 

e(z + r ) e(z + n e(z + 

(2.9) D„i,,lne(z) = „+<„5<^±|||^, 

for any z G and any distinct points a,b,c G TZg] here the scalars Pi,qi, i = 1,2 depend 
on the points a,b,c and are given by 

(2.10) p,(a,6,c) = -Aln^|^, 

(2.11) ^^("'^'^^^eSioD''^^^'^^"^' 

(2.12) qi{a,b) = DaDblnemf'' 

(2.13) 92(0,6) 



a • 



Dae[6]{0)DbQm0) 



may 

where (5 is a non-singular odd characteristics. 
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2.2. Real Riemann surfaces. A Riemann surface Rg is called real if it admits an an- 
tiholomorpliic involution, denoted by r. An anti-liolomorpliic involution r : TZg — )■ TZg 
satisfies = id and acts on the local parameter as the complex conjugation. The con- 
nected components of the set of fixed points of the anti-involution r are called real ovals 
of T. We denote by TZg{M.) the set of fixed points. According to Harnack's inequality 
|26) . the number x of real ovals of a real Riemann surface of genus g cannot exceed g + 1: 
< X ^ 5 + 1- Curves with the maximal number x = 9 + 1 of real ovals are called M-curves. 

The complement TZg\TZg{M.) has either one or two connected components. The curve TZg 
is called a dividing curve if TZg \ TZg{M.) has two components, and TZg is called non-dividing 
if TZg \ TZg{M.) is connected (notice that an M-curve is always a dividing curve). In this 
paper we only consider hyperelliptic curves which are discussed in detail section 4. 

Let {A,B) be a basis of the homology group Hi{TZg). According to Proposition 2.2 in 
Vinnikov's paper |32j (see also [25]), there exists a canonical homology basis such that 



(2.14) 



tA 
tB 







where I„ is the g x g unit matrix, and H is a block diagonal g x g matrix defined as follows: 



1) ifTZgiM.) / 



/O 1 
1 



1 

1 



if TZg is dividing, 



0/ 



/I 



V 

rank(]HI) = g + 1 — x both cases. 



if TZg is non-dividing; 



0/ 



2) HTZgil 



(i.e. the curve does not have real ovals), then 



/O 1 
1 



V 



1 

1 0/ 



or 



/O 1 
1 



V 

rank(BI) = g if g is even, rank(]HI) = g — 1 if 5 is odd. 



1 

1 



0/ 
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In what follows we choose a canonical homology basis in HiilZg) satisfying (2.14) and 
take a,b G TZg such that ra = b. Denote by ^ a contour connecting the points a and 
b which does not intersect the canonical homology basis. Then the action of r on the 
generators (A,B,£) of the relative homology group Hi{TZg, {a,b}) is given by (see [27] for 
more details) 



(2.15) 




Ig 

M -Ig 

N* 




for some N G Z^. In the case where ra = a and rb = b, the action of r on the generators 
{A,B,i) of the relative homology group Hi{TZg, {a,b}) reads (see |27| for more details) 

(2.16) I tB 




Ti J \ N* 

where the vectors N, M G are related by 

(2.17) 2N + ]HIM = 0. 

Now let us study the action of r on Abelian differentials and the action of the complex 
conjugation on the theta function with zero characteristics. Denote by r* the action of r 



lifted to the space of differentials: T*uj{p) = io{Tp) for any p G IZg. By (2.14) the ^-cycles 
of the homology basis are invariant under r. Due to the normalization conditions (2.1), 
this leads to the following action of r on the normalized holomorphic differentials: 



(2.18) T*ujj = -ujj. 

Let a,& G TZg and denote by {A,B^Sh} the generators of the homology group Hi{TZg \ 
{a, 6}) of the punctured Riemann surface lZg\{a, 6}, where Sh is a positively oriented small 
contour around b such that S^o i = 1. It was proved in [27] that in the case where ra = 6, 
the ^-cycles in the homology group Hi{lZg \ {a, 6}) are stable under r. Therefore, by the 
uniqueness of the normalized differential of the third kind Q-b-a which has residue 1 at 6 
and residue —1 at a, we get 



(2.19) T*nb-a = -^b- 



a- 



In the case where ra = a and rb = b, Proposition A. 2 in [27] shows that the action of r 
on the ^-cycles in the homology group Hi{TZg \ {a, b}) is given by 

(2.20) tA = A- MSb, 



where M is defined in (2.16). Therefore, by the uniqueness of the differential ftb-a, we 
deduce that 

(2.21) T*nb^a = ^b-a + M'OJ, 



where to denotes the vector of normalized holomorphic differentials. From (2.14) and (2.18) 
we obtain the following reality property for the Riemann matrix B: 

(2.22) l = ]B-2i7r]H. 



Moreover, according to Proposition 2.3 in [32j, for any z G C^, relation (2.22) implies 



(2.23) e(z) = Ke(z - i7rdiag(M)), 
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where diag(]H) denotes the vector of diagonal elements of the matrix H, and k is a root of 
unity which depends on the matrix H (knowledge of the exact value of n is not needed for 
our purpose). 



2.3. Action of r on the Jacobian and the theta divisor of real Riemann surfaces. 

In this subsection, we review known results about the theta divisor of real Riemann surfaces 



(see |32| 116)). Let us choose a canonical homology basis satisfying (2.14) and consider the 
Jacobian J := J{TZg) of a real Riemann surface TZg. 

The anti-holomorphic involution r on TZg gives rise to an anti-holomorphic involution 
on the Jacobian: if 2? := 2?2 — with Pi and D2 positive divisors on TZg (recall that 
a positive divisor is defined by a finite formal sum of points Y^- rii a, with Oj G TZg and 
Ui G N) then rP is the class of the point {J^^^ uj) = {J-p^ t*Ljj) in the Jacobian. Therefore, 



by ( |2.18 ) r lifts to the anti-holomorphic involution on J, denoted also by r, given by 



(2.24) rC = -C, 
for any C G J. 

Now consider the following two subsets of the Jacobian 

(2.25) 5i = {CG J; C + rC = ivrdiag(H)}, 

(2.26) 52 = {CG ^; C--rC = ivrdiag(]H)}, 



where the matrix HI was introduced in ( 2.14^ . Below we study their intersections Si R (0) 



and 5*2 D (0) with the theta divisor (0), the set of zeros of the theta function. Let us 
introduce the notation: the vectors a, i = 1, ... ,5 with components e^fc = 6ik, Bj = Be,. 

It is a straightforward computation to prove that the set Si is the disjoint union of the 
tori defined by 

(2.27) 

r„ = {C G J; C = ivr (diag(H)/2 + vi Cr+i + ... + Vg-r eg) + /3i Re(Bi) + . . . + ^g Re{Mg) , 

/3i,...,/3, GM/2Z, /?,+!,..., gM/Z}, 

where v = {vi, . . . ,Vg-r) G (Z/2Z)3~'' and r is the rank of the matrix H. Therefore, the 
description of the set Si D (0) reduces to the study of the sets D (0). In the case where 
TZg{M) 7^ and TZg is non-dividing, one can see that for all v the torus contains a 
half-period corresponding to an odd half-integer characteristics, which yields D (0) 7^ 0. 
The same holds for all u 7^ in the case where the curve is dividing or does not have 
real ovals. The following result proved in [32] provides a complete description of the sets 
Ty n (0) in the case where the curve admits real ovals: 

Proposition 2.1. IfTZg{M) 7^ 0, then n (0) = ^ if and only if the curve is dividing and 
V = 0. 

In other words, among all curves which admit real ovals, the only torus Ty which does not 
intersect the theta divisor is the torus Tq corresponding to dividing curves. This torus is 
given by 
(2.28) 

To = {C G J; C = /3i Re(Bi) + . . . + /3<, Re{Bg), A, . . . , /3, G M/2Z , /3,+i, . . . , /^^ G R/Z}. 
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Analogously, it can be checked that the set S2 is the disjoint union of the tori defined 

by 

(2.29) n = {C € J ; C = 2m {aiei + . . . + ag eg) + {vi/2) B,+i + . . . + {vg^r/2) Mg, 

ai,...,ageR/Z}, 

where v = (fi, . . . ,Vg-r) € (Z/2Z)^~^ and r is the rank of the matrix H. Description of 
the sets D (0) in the case where the curve admits real ovals was given in |16]: 

Proposition 2.2. IfTZgCR) ^ 0, then n (0) =% if and only if the curve is an M-curve 
and V = 0. 

3. Algebro-geometric solutions of the Camassa-Holm equation 

In this section we will use Fay's identities to construct solutions to the CH equation 
on hyperelliptic surfaces. For the resulting formulae we establish conditions under which 
we obtain real and smooth solutions. In what follows TZg denotes a hyperelliptic curve of 
genus g > 0, written as 

2g+2 

(3.1) i,^=ll{X-X,), 

where the branch points A, G C satisfy the relations Aj 7^ Xj for i ^ j- We denote by a 
the hyperelliptic involution defined by a{X,fi) = (A, —fJ.)- Note that the CH equation can 
be expressed in the following simple form, 

(3.2) mt + urux + 2mux = 0, 
where we put m := u — Uxx + k. 

3.1. Identities between theta functions. In our approach to construct algebro-geometric 
solutions of the CH equation we use the corollaries (2.8) and (2.9) of Fay's identity. 

Proposition 3.1. Let a,b £ TZg such that cr{a) = b and let e £ TZg be a ramification 
point, namely, e = {Xj,0) for some j £ {1, . . . ,2g + 2} . Denote by gi and g2 the following 
functions of the variable z G C^; 

where r = J uj and u is the vector of normalized holomorphic differentials. Then the two 
following identities hold: 

(3.4) DbDeln^ = -^Dblngig2, 

92 9192 

(3.5) DhDeln{gig2) = ^ ( Dfcln— -2pi] - 2925152- 

P2 gm \ g2 J 

Here we used the notation: 

(3.6) p2 = P2{b,e,a), pi = pi{e,b, a), q2 = q2{b,e), 

where the scalars q2{-,-) and pi{., .), i = 1,2, are defined in {2A^ and \2.10 ), (2.11 ); 
(respectively De) denotes the directional derivative along the vector V;, (respectively Vg) 
defined in (2.6). 
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Proof. Under the changes of variables {a,b,c) — t- (6, e, a) and z — )■ z — r/2, identity (2.8) 
becomes 



(3.7) 



D, In 



91 



- Pi + 



P2 



92 9192 

where we used the notation pi = pi{b,e,a) for i = 1,2. Here we used the fact that 
(jj = J u! = r/2, according to the action of a* (the action of a hfted to the space of 
one-forms) on the normahzed holomorphic differentials ajj: 



(3.8) 



-UJj, 



J 



,9- 



Applying the differential operator Di, to equation (3.7) one gets: 

Db{9i92) 



DbDeln 



91 

92 



-P2 



{9i92f 



P2 

9192 



Dh ln(/i52, 



which proves (3.4). To prove (3.5), consider the change of variables {a,b,c) 
(2.8), which leads to 

n 1 ~ ^~ 0(z + r) 

Dblngi = pi +P2 92 

Changing z to — z in the last equality, one gets 

Di,\ng2 = -pi -P29i 



(e, b, a) in 



e(z + |)- 

9(z - r) 
Gfz-^) 



From these two identities, it can be deduced that 

e(z + r) 



(3.9) 
(3.10) 

Moreover, since 



e(z + |) 
e(z - r) 



{P292) ^ {DblTLQi -pi), 

- {p2 9iy^ {Db 11152 +pi] 



DbDeHgig2) = DbDelnQ(z + - ) + DbDelnOlz 



using (2.9) one gets 
(3.11) 



DbDe ln{gig2) 



q2 0(z + r) ^q2 6(z - r) 



5iG(z+§) g2e(z-|) 



2 DbDe Ine(z) 



2 92(^152, 



which by (3.9) and (3.10) leads to B.bh. 



□ 



With identities (3.4) and (3.5) we are now able to construct theta-functional solutions of 
the CH equation: 

Theorem 3.1. Let a,b ^ TZg such that cr{a) = b, and let e £ TZg be a ramification point. 
Denote by I an oriented contour between a and b which contains the point e. Assume that 
i does not cross cycles of the canonical homology basis. Choose arbitrary constants d € O 
and k, G C, and put 

(3.12) ai = pi{b,e,a), 02 = 2pi{e,b,a) + k, 



where the function pi is defined in (2.10). Let y{x,t) be an implicit function of the variables 
X, t G M defined by 

e(z-d + |) 



(3.13) 



X + aiy + a2t + C = 



efz-d 



)' 
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where r = J^oj. Here the vector Z is given by 

(3.14) Z{x,t)=Vey{x,t) + \bt, 

where the vectors Vg and V;, are defined in ^2^. Then the following function of the 
variables x and t is solution of the CH equation: 



(3.15) 



^ ^ e(z-d + |) 

U{X, t) = Db In ; ^ - 02- 



G(Z-d-0 

Here denotes the directional derivative along the vector V;,. 

Note that the function y{x,t) is the same function as introduced in |[3|. 

Proof. Let P,6 £ C and ai, 02 £ C he arbitrary constants. Let us look for solutions u of 
CH having the form 

(3.16) uix,t) = (3Dhln®^^^^^+5 = l3Dbln^ + 6, 



where Z(x,t) is defined in (3.14), and the functions gi,g2 were introduced in (3.3) with 
z = Z(x,t) — d. By (3.13), the derivative with respect to the variable x of the implicit 



function y{x,t) is given by 



DAn 



91 

92 



ai 



With ai = pi{b,e,a), relation (3.7) implies 

(3.17) = 
Analogously it can be checked that 

(3.18) y^ = -y^ 



_ 9192 

P2 

u 5 



0.2 



Now let us express the function m{x, t) = u — Uxx + k introduced in (3.2 ) in terms of the 
functions gi and g2 of (3.3). By (3.4) and (3.17), the first derivative of the function u 
( 3.16| ) with respect to the variable x is given by 

(3.19) Ux = -pDbHgig2). 



By (3.5) and (3.17) we obtain for the second derivative of u with respect to x: 
(3.20) = (^1)6 In ^ - 2p^ - 2(3 p2 {gm? ; 



here we used the identity q2 = —P2P2 relating the scalars ^"2,^2 and p2 defined in (3.6). 
Therefore, with (3.16) and (3.20), the function m reads 

(3.21) m(x, t) = 6 + k + 2l3pi + 2pp2 {gmf- 

Taking the derivative of m with respect to x, and the derivative of m with respect to t, 
one gets respectively: 

(3.22) rnxix,t) = 4.(3 p2{gig2f Vx De\n{gig2), 

(3.23) mt(x,t) = 4/3 p2 (5152)^ (yt L>e ln(5i52) + ln((7i52)) • 
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Therefore, substituting the functions (3.16), (3.19), (3.22) and (3.23) in the left-hand side 



of the CH equation (3.2) we obtain 



2p2Dehi{gig2)yx (gm)^ n ^1 - -^^ + ^ + 02 

-Dbln{gig2) {6 + k + 2/3 p, + 2p2 (gmf W - I)) = 0. 
This equality holds for /3 = 1, 5 = — 2pi — k and 02 = — <5, which completes the proof. □ 



3.2. Real-valued solutions and smoothness conditions. In this subsection, we iden- 
tify real-valued and smooth solutions of the CH equation among the solutions given in 
Proposition 3.2. Let us first recall that hyperelliptic M-curves of genus g can be given by 
the equation 



(3.24) 



2g+2 



A. 



i=l 



where the branch points Aj are real and satisfy Aj 7^ Aj i ^ j. On such a curve, we 
can define two anti-holomorphic involutions ri and T2, given respectively by ti{X, ^) = 
{X,Ji) and T2(A,^) = {X,—Jl). Let us show that the curve (3.24) is an M-curve with 
respect to both anti-involutions ri and T2- In the case where Aj G M satisfy Ai < ... < 
A2g+2) it can be seen that projections of real ovals of ri on the A-plane coincide with the 
intervals [A2g+2) Ai], [A2, A3], . . . , [A2g, A2g+i], whereas projections of real ovals of T2 on the 



A-plane coincide with the intervals [Ai, A2], . . • , [A2g+i, A2g+2]- Hence the curve (3.24) has 
the maximal number g + 1 oi real ovals with respect to both anti-involutions ti and T2. 
Now assume that TZg is a real hyperelliptic curve which admits real ovals with respect to 



an anti-holomorphic involution r. Let us choose a homology basis satisfying (2.14). Recall 



that TZg{M.) denotes the set of fixed points of the anti-holomorphic involution r. 

The following propositions provide reality and smoothness conditions for the solutions 



u{x, t) (3.15 ) in the case where the points a and b are stable under r, and in the case where 
TO = b. It is proved that, for fixed to ^ the function ti(x,to) is smooth with respect to 
the real variable x when a and b are stable under r. In the case where ra = b, the function 
u{x, to) is either smooth, or it has cusp-like singularities. 

Proposition 3.2. Assume that TZg is a hyperelliptic M-curve of genus g, and denote by 
e G TZgiM.) one of its ramification points. Let a,b £ TZg{M.) such that cr(a) = b. For any 
c £ {a,b,e}, choose a local parameter kc such that kc{Tp) = kc{p) for any point p in a 
neighborhood of c. Denote by i an oriented contour between a and b containing point e 
which does not intersect cycles of the canonical homology basis. Choose I such that the 

Choose 



closed path ri 

C G C in (3.13^ such that /m(C) 



is homologous to zero in Hi{lZg). Take d G iM^ and k G 



arg 



In 



Then solutions u{x, t) of the 
CH equation given in (3.15) are real-valued, and for fixed to G M, the function n(x, to) is 



e(d+r/2) 
0(d-r/2) 



smooth with respect to the real variable x. 



Proof. Let us check that under the conditions of the proposition, the function u{x, t) (3.15 ) 
is real-valued. Let us fix y, t G M. First of all, invariance with respect to the anti-involution 
T of the points e and b implies 



(3.25) 
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where the vector Z is defined in (3.14). In fact, using the expansion (2.6) of the normahzed 
holomorphic differentials loj near c G {e, 6}, one gets 

for any point p in a neighborhood of c. Then by ( |2.18[ ), the vectors Vg and Vf, appearing 
in the vector Z are purely imaginary, which leads to (|3.25 ). Moreover, since the closed 
contour ri — i is homologous to zero in Hi(TZg), from (2.18) one gets 

(3.26) r=-r. 

For arbitrary points 01,02,03 € TZg, using the representation of the differential Qjj^—^jj^ 
terms of multi-dimensional theta functions (see, for instance, [^), one gets: 



m 



(3.27) 



^as—ai {p) 



dkaiip) 



p=a2 



pi{ai,a2, as). 



We deduce that the scalars pi{b,e,a) and pi{e,b,a) appearing respectively in ai and 02 
(see ( |3.12 )) satisfy 



(3.28) 



^a-b{p) 



dke{p) 



pi{b,e,a), 



^a-e{p) 



p=b 



pi{e,b,a). 



p=e dkb{p) 

Therefore, since the points a,b,e are stable under r, from (2.21) it can be deduced that 
Pi{b, e, a) and pi{e, b, a) are real, which involves 

(3.29) oi = «!, 02 = a2- 

Let y,t GM and denote by h the function 

e(z-d + |) 



(3.30) 



h{y,t) 



e z-d 



By (3.13), X is a real- valued function of the real variables y and t if the function h is real 
and has a constant sign, and if we choose Im((^) = arg{ln(/i(0, 0))}. From (2.23), (3.25) 
and (3.26) we deduce that 

(3.31) h{y,t) 



e(Z + d+|+i7rdiag(]H)) 
e(Z + d- § + i7rdiag(IHI)) 
Let us choose a vector d G such that 

d = - d - ivr diag(H) + 2mT + BL 



for some vectors T, L G Z^. Reality of the vector d -|- d together with (2.22) implies 
(3.32) d = -Re(]B)L -Fid/ 



for some d/ S and the relation 2 T -|- HL = diag(]HI) for T and L. For this choice of the 

e(z-d + |) 



vector d, (3.31) becomes 



h{y,t) 



efz-d 



exp{- (r,L)}, 



where we used the quasi-periodicity property (2.4) of the theta function. Therefore, the 
function h is real if L = 0, that is d E iM^. Now let us check that h has a constant 

the functions 



d ± I G iM^, by Proposition 



2.2 



0, that is d E iRf 
sign with respect to y,t £ M. Since Z 

e(z 

curve, i.e., if all branch points in (3.1 ) are real. Hence /i is a real continuous non vanishing 



d lb I ) of the real variables y and t do not vanish if the hyperelliptic curve is an M- 
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function with respect to the real variables y and t, which means it has a constant sign. 



Therefore, x is a real- valued function of y and t if the constant C in (3.13) is chosen such 
that Im(C) = arg{ln(/i(0, 0))}. It is straightforward to see that the solution u (3.15) is a 



real- valued function of the real variables y and t, and then is real- valued with respect to 
the real variables x and t. 

Now fix to S M and let us study smoothness conditions for the function u{x, to) with 



respect to the variable x. First let us check that the solution u (3.15) is a smooth function 



of the real variable y if it does not have singularities. Since the theta function is entire, 
singularities of the solution u are located at the zeros of its denominator. As we saw in the 
previous paragraph, if the curve is an M-curve and d G iM^, the functions 0(Z — d it |) 
and B(Z — d) do not vanish. In this case, the function u is smooth with respect to the 
real variable y. Now let us prove that u is smooth with respect to the real variable x. By 
( 3.13| ), the function x{y) is smooth. Moreover, it can be seen from (3.13) and (3.7) that 

e(z-d)2 



(3.33) 



P2 



9i( 



■ P2 



e z 



e z-d 



-d)52(Z-d) 

Since the functions 0(Z — d it |) and B(Z — d) do not vanish, we deduce that x{y) is a 
strictly monotonic real function, and thus the inverse function y{x) has the same property. 
Therefore, the function u{x,tQ) = u{y{x)) is a smooth real-valued function with respect to 
the real variable x. □ 

Now let us study real-valuedness and smoothness of the solutions in the case where ra = b. 

Proposition 3.3. Assume that TZg is a hyperelliptic M-curve of genus g, and denote by 
e G TZg{M) one of its ramification points. Let a,b € TZg such that a{a) = b and assume that 
TO = b. Choose local parameters such that kb^rp) = ka{p) for any point p in a neighborhood 
of a, and ke{Tp) = — ke{p) for any p lying in a neighborhood of e. Denote by i an oriented 
contour between a and b containing point e, which does not intersect cycles of the canonical 
homology basis. Assume that N = 2L for some L S Z^, where 'N € is d efined in (2.15). 
Take A; € M and define d = d/j+ y N for some d/j G R^. Choose £ C in (3.13) such that 

Im{C) = arg 



|ln (^ 0(^-1-/2) ) }• ^^6"- solutions u (3.15) of the CH equation are real-valued. 



Moreover, for fixed tQ G M, the function u{x, to) is smooth with respect to the real variable 
X in the case where N = 0, otherwise it has an infinite number of singularities of the type 

O [{x — a;o)2"+i J for some n G N \ {0} and xo G M, i.e., cusps. 

Proof. Analogously to the case where a and b are stable under r, let us prove that solutions 
u ( 3.15[ ) are real-valued. Fix y, t G M. First let us check that the vector Z (3.14) satisfies: 

(3.34) Z = Z. 

From (2.6) and (2.18) one gets Vq = — and Vg = Vg. Moreover, the vectors Vq and 
Vfc satisfy -|- = because of (3.8 ); thus we have = and Vg = Ve which proves 



( 3.34[ ). By (2.19) and (3.28) it can be deduced that a\ and 02 (3.12) satisfy 
(3.35) oi = ai, 02 = «2- 



Moreover, from (2.18) and (2.15) one gets: 
(3.36) i 



r - 2i7rN, 



where N G is defined in (2.15). Let us check that x (3.13) is a real-valued function of the 



real variables y and t. By (3.35 ), this holds if the function h (3.30 ) is real and has a constant 
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sign with respect to the real variables y and t, and if we choose Im(C) = arg{ln(/i(0, 0))}. 
By (|2.23D, KM and (|3.36[) it follows that 



h{y,t) 



G(Z-d + §+p) 



(3.37) 
where p 

which is, since d — d and p are purely imaginary, equivalent to d = d + p + 2i7rT, for 



e(z-d-| + p)' 

IttN — i7rdiag(BI). Let us choose the vector d G such that 
d = d + p (mod 2i7rZ^^ + : 



some T G Z^. Here we used the action (2.22) of the complex conjugation on the Riemann 
matrix B, and the fact that B has a negative definite real part. Hence, the vector d can 
be written as 



(3.38) 



ITT 

d = di^+-(N + diag(]H) 



2T) 



for some d^ G and T G Z^. For this choice of the vector d, by (3.37) the function h 



is a real-valued function of the real variables y and t. Now let us study in which cases the 
function h has a constant sign. The sign of the function h is constant with respect to y 
and t if the functions 0(Z 

Z ^ 

0(Z — d lb |) do not vanish if the hyperelliptic curve is dividing (in this case diag(IHI) = 0) 
and if the arguments Z 



do not vanish. By (|3.34|), (|3.36|) and (|3.38), the vectors 

the functions 



-dib I belong to the set Si introduced in (2.25 ). Hence by Proposition 
di"" 



2.1 



' d + ^ is real if T 



d z fc I in the theta function are real (modulo 2i7rZ3). The vector 



in (3.38). With this choice of the vector T, the imaginary part of 



the vector Z — d — | equals — ivrlN . Therefore, the vector Z — d — | is real modulo 2i7rZ^ if all 



components of the vector N are even. To summarize, the function h{y, t) defined in (3.30) 
is a real- valued function with constant sign if the hyperelliptic curve is dividing (i.e., all 
ramification points are stable under r, since the ramification point e is stable under r and 
since dividing curves have either only real branch points, or pairwise conjugate ones), if 
T = and N = 2L for some L G Z^, where vector N G Z^ is defined in (2.15 ). Analogously 



to the proof of Proposition 3.2 we conclude that x is a real- valued continuous function of 



the real variables y and t, and thus solutions u{x,t) (3.15) are real- valued functions of the 
real variables x and t. 

Now let us study smoothness conditions for fixed to G M. Notice that the function u{y) 



(3.15) is a smooth function of the real variable y since the denominator does not vanish. 



as we have seen before. Put z = Z — d. Let us consider the function Xy{y) given in (3.33) 
in both cases: N = and N 7^ 0. 

- If N = 0, the function Xy{y) does not vanish, since in this case z G which implies 
that the function 0(z) does not vanish. Hence, analogously to the case where a and b are 
stable under r, for fixed to S '^1 the function ti(x,to) is smooth with respect to the real 
variable x. 

- If N 7^ 0, the function 0(z) vanishes when z belongs to the theta divisor. Fix Xo,to S M 
and denote by zq and yo the corresponding values of z and y. Assume that zq is a zero of 



the theta function of order n > 1. Then by (3.33), the function Xy{y) has a zero at yo of 
order 2n. It follows that function x{y) — x{iiq) has a zero of order 2n -|- 1 at yoj and then 



(3.39) 



y{x) - yo = O ((x - xo) 
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On the other hand, it can be seen from (3.4) that 

e(z) 



(3.40) 



Uy{y) =P2 



e(z + |)e(z-|) 



[e(z)v(z)-2L»be(z) 



where ijj{z) = Dhln (0(z + |) 0(z — |)). Identity (3.40) imphes that function Uy{y) has a 
zero at yo of order 2n — 1, namely, 

(3.41) u{y)-u{yo) = 0{{y-yof''). 



Finally with (3.39) and (3.41 ), the function u(x, to) has an infinite number of singularities 
of the type O (^{x — xq) j , i.e., cusps. □ 

4. Numerical study of algebro-geometric solutions to the Camassa-Holm 

equation 



In this section we will numerically study concrete examples for the CH solutions (3.15). 



As shown in the previous sections, real and bounded solutions are obtained on hyperelliptic 
M-curves, i.e., curves of the form 

2g+2 



n - 



where g is the genus of the Riemann surface, and where we have for the branch points 
Aj E M the relations Aj ^ Xj for i ^ j- 



For the numerical evaluation of the CH solutions (3.15) we use the code presented in 



|201l2T| for real hyperelliptic Riemann surfaces. The reader is referred to these publications 
for details. The basic idea is to introduce a convenient homology basis on the related 
surfaces, see Fig. [T] It is related to the basis used in the previous sections by the simple 




:g+2 



Figure 1. Homology basis on real hyperelliptic M-curves, contours on sheet 
1 are solid, contours on sheet 2 are dashed. 

relation A — )■ —B, B ^ A. This choice of the homology basis has the advantage that 
the limit in which branch points encircled by the same ^-cycle collide can be treated 
essentially numerically. Below we will consider examples where the distance of such a pair 
of branch points is of the order of machine precision (10-^"^) and thus numerically zero. 
This limit is interesting since the ;B-periods diverge, which implies that the corresponding 
theta functions reduce to elementary functions. The CH solutions (3.15) reduce in this 



case to solitons or cuspons. Since we want to study also this limit numerically, we use the 
homology basis of Fig. [T] and not the one of the previous sections. 

The sheets are identified at the point a by the sign of the root picked by Matlab. We 
denote a point in the first sheet with projection A into the complex plane by A^^^ and 
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a point in the second sheet with the same projection by A^^^. The theta functions are 
in general approximated numerically by a truncated series as explained in [20] and [Sj. A 
vector of holomorphic differentials for these surfaces is given by (1, A, ... , A^~^)* dX/fi. The 
periods of the surface are computed as integrals between branch points of these differentials 
as detailed in \21\ . The Abel map of the point a (and analogously for b) is computed in a 
similar way, see |23, as the integral between a and the branch point with minimal distance 
to a. It is well known (see for instance |7|) that the Abel map between two branch points 
is a half period. 

To control the accuracy of the numerical solutions, we use essentially two approaches. 
First we check the theta identity (2.8), which is the underlying reason for the studied 
functions being solutions to CH. Since this identity is not built into the code, it provides 
a strong test. This check for various combinations of the points a,b and e ensures that 
the theta functions are computed with sufficient precision, and that the quantities ai and 
02 in (3.13) are known with the wanted precision (we always use machine precision here). 
In addition, the smooth solutions are computed on Chebyshev collocation points (see, 
for instance, ^31j) for x and t. This can be used to approximate the computed solution 
via Chebyshev polynomials, a so-called spectral method having in practice exponential 
convergence for smooth functions. Since the derivatives of the Chebyshev polynomials can 
be expressed linearly in terms of Chebyshev polynomials, a derivative acts on the space 
of polynomials via a so called differentiation matrix. With these standard Chebyshev 
differentiation matrices (see |3T]), the solution can be numerically differentiated. The 
computed derivatives allow to check with which numerical precision the partial differential 
equation (PDE) is satisfied by a numerical solution. With these two independent tests, 
we ensure that the shown solutions are correct to much better than plotting accuracy (the 
code reports a warning if the above tests are not satisfied to better than 10~^). We do not 
use the expansion in terms of Chebyshev polynomials for the cusped solutions since the 
convergence is slow for functions with cusps. 

We first consider smooth solutions u (3.15) in the case ra = a, rb = b. To obtain non- 
trivial solutions in the solitonic limit, we use a vector d corresponding to the characteristics 
5 [o " of examples. To plot a solution u in dependence of x and t, we compute it on 

a numerical grid for y and t to obtain x{y,t) defined in (3.13) and u{y,t) given by (3.15). 
These are then used to obtain a plot of u{x, t) without having to solve the implicit relation 
(3.13). In all examples we have k = 1. The solutions for ra = b, e.g., for points a and b on 
the cuts encircled by the ^-periods in Fig. [T] look very similar and are therefore not shown 
here. 

Solutions on elliptic surfaces describe travelling waves and will not be discussed here. In 
genus 2 we obtain CH solutions of the form shown in Fig. [2] The typical soliton collision 
known from the KdV equation is also present here, the unchanged shape of the solitons 
after the collision, but an asymptotic change of phase. 

In genus 6 the CH solutions have the form shown in Fig. [3] In the solitonic limit one 
can recognize a 6 soliton event. 

The reality properties of the quantities entering the solution ( |3.15[ ) depend on the choice 
of the homology basis. For instance, in the homology basis of Fig. [Ijand for a and b stable 
under r, the Abel map r up to a vector proportional to in and the vectors V^, Vg are real, 
whereas these quantities are purely imaginary in the homology basis used in the previous 
sections. Thus the easiest way to obtain cusped solutions is in this case to put e = A2 and 
to choose d corresponding to the characteristics 5 [} i]*- It can be easily checked that 
the theta functions B(Z — d it r/2) cannot vanish since the argument is real, whereas the 
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Figure 2. Solution (3.15) to the CH equation on a hyperelliptic curve of 
genus 2 with branch points —3, —2, 0, e, 2, 2 + e and a = (—4)*^"'^^, b = (— 4)*^^^ 
and e = (—3,0) for e = 1 on the left and e = 10~^^, the almost solitonic 
limit, on the right. 





Figure 3. Solution (3.15) to the CH equation on a hyperelliptic curve of 
genus 6 with branch points —7, —6, —5, —5 + e, —3, —3 + e, — 1, — 1 + e, 1, 1 + 
e, 3, 3 + e, 5, 5 + e ari(i a = (-8)(i), b = (-8)(2) and e = (-7, 0) /or e = 1 on 
the left and e = 10~^^, the almost solitonic limit, on the right. 



0(Z — d) will have zeroes since the argument is complex. This implies that the derivative 
Ux in (3.19) diverges which corresponds to cups for the solution u. Peakons do not appear 
in such a limit of theta-functional solutions to CH and are thus not discussed here. To 
obtain them one would have to glue solutions in the solitonic limit on finite intervals to 
obtain a continuous solution that is piecewise C^. 

We show the cusped solutions always in a comoving frame x' = x + vt to allow a better 
visualization of the solutions. In genus 2 we obtain cusped CH solutions of the form shown 
in Fig. [4j where also cusped solitons can be seen in the degenerate situation. Obviously 
the collision between cuspons is analogous to soliton collisions. 

In genus 6 the CH solutions have the form shown in Fig. [5] In the solitonic limit one 
can recognize a 6-cuspon. 
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Figure 4. Cusped solution (3.15) to the CH equation on a hyperelliptic 
curve of genus 2 with branch points — 3, — 2, 0, e, 2, 2 + e and a = (—4)^^^, 
b = (—4)*^^^ and e = (—2,0) for e = 1 on the left and e = 10~^^, the almost 
solitonic limit, on the right. 



Figure 5. Solution (3.15) to the CH equation on a hyperelliptic curve of 
genus 6 with branch points —7, —6, —5, —5 + e, —3, —3 + e, — 1, — 1 + e, 1, 1 + 
e, 3, 3 + e, 5, 5 + e and a = (-8)(i), b = (-8)(2) and e = (-6, 0) /or e = 1 on 
the left and e = 10""^^, the almost solitonic limit, on the right. 



5. Conclusion 

In this paper we have shown at the example of the CH equation that Mumford's program 
to construct algebro-geometric solutions to integrable PDEs can be also applied to non-local 
(here in x) equations. For the studied case the solutions in terms of multi-dimensional theta 
functions do not depend directly on the physical coordinates x and t, but via an implicit 
function. One consequence of this non-locality is the existence of non-smooth solitons. A 
numerical study of smooth and non-smooth solutions was presented. 

A further example in this context would be the equation from the Dym-hierarchy for 
which theta-functional solutions were studied in [2], which will be treated elsewhere with 
Mumford's approach. It is an interesting question whether the 2+1 dimensional generaliza- 
tion of the CH equation |17| 114]. for which algebro-geometric solutions are so far unknown, 
can be also treated with these methods. 
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